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Abstract. Observations of solar prominences reveal a complex, dynamic flow 
field within them. The flow field within quiescent prominences is characterized 
by long "threads" and dark "bubbles" that fall and rise (respectively) in a thin 
sheet. The flow field in active prominences display more helical motions that 
travel along the axis of the prominence. We explore the possible dynamics of 
both of these with the aid of 2.5D MHD simulations. Our model, compressible 
plasma possesses density and temperature gradients and resides in magnetic 
field configurations that mime those of a solar prominence. We present results 
of various configurations and discuss the nonlinear behavior of the resulting 
dynamics. 



1 Introduction 



Recent observations of solar prominences have revealed a complex, dynamic flow 
field within them. The flow field within quiescent prominences is characterized 
by long "threads" and dark "bubbles" that fall and rise (respectively) in a thin 
sheet. The flow field in active prominences display more helical motions that 
travel along the axis of the prominence. We investigate these behaviors with the 
aid of axisymmetric MHD simulations. 

Flows within prominences may be influenced by double diffusive instabili- 
ties, whic h have been ext ensively studied for the role they may play in the solar 
dynamo Hughes fc Weiss | ([1995). In this case the magnetic fields provide buoy- 
ancy to a rising a fluid parcel through the magnetic pressure even if the thermal 

background is stably stratified. This can occue when dB d J p < 0, where B is the 
horizontal magnetic field strength, p is the density and z is the height. In other 
words, when the Alfven speed in a fluid layer decreases with height, magnetic 
"fingering" instabilities are possible. 

It is quite likely that such conditions can be met in prominences. If promi- 
nences are plasma trapped about the solar photosphere by a helical field whose 
axis is aligned to the horizontal, then the region below the axis of the flux tube 
will indeed be a location where the field decreases with height. If this corre- 
sponds to a region with relatively large density scale heights, fingering instabili- 
ties are possible. Here we concoct a simple axisymmetric model of a prominence 
to investigate the feasibility and possible signatures of any such dynamics. We 
present our model in the next section, display our numerical results in the next 
and close with a discussion and conclusion. 
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Figure 1. The detailed properties of one buoyant plume at mid-rise shows 
that the buoyancy due to the effect of magnetic pressure rather than temper- 
ature. The result is a cool, low-density object that would appear as a dark 
feature Hinode images. The radial (vertical) velocity u r and entropy S are 
also shown. 



2 Model 



Partial differential equations describing compressible MHD are solved in an ax- 
isym metric cylindrical geometry, u si ng a numerica l code developed for this pur- 
pose (|Hurlburt fc Rucklidge||2000|). |Botha et af] fl2008D extend this model to 



include azimuthal flows and fields. Here we modify the initial state and back- 
ground stratification to mimic a horizontal, helical flux tube and consider the 
dynamics of buoyancy instabilities. 

The computational domain is an axisymmetric cylinder of inner, outer ra- 
dius Ti,r2, situated in the (r,z) plane. At the radial boundaries boundary the 
magnetic field is imposed to be consistent with its initial condition; the tempera- 
ture gradient is chosen to be zero and the horizontal boundaries are impenetrable 
and stress free. The inner and outer walls are perfect conductors and field lines 
are tied on the surfaces z = 0, 1. The plasma is triggered by starting the quiet, 
hydrostatic equilibrium perturbed by a hot blob in pressure equilibrium near 
the outer boundary. The initial temperature and density profiles in the radial 
(r) direction are T = T , p = p e^ 2/2 , = B (r/r 2 ) 2 sin (9, and B z = B cos9. 
Here, the subscript defining the quantity at the outer edge of the cylinder 
(r = IT^); 9 represents the initial angle of the field twist along the outer wall; 
and the gravitational force constant g~ x — HTqT2 where H is the pressure scale 
height at r = T2. 

The physical quantities are dimensionless, with the length scaled by the 
length of cylinder, velocity scaled by the sound speed in it and temperature, 
magnetic field, density, and pressure all scaled to their initial values at r = 0. 



3 Numerical results 

With a purely azimuthal field, the plumes form and rise with classic plume 
profile as seen in Figure [Tj The buoyancy driving from the initial thermal pulse 
is replaced by the magnetic pressure and produces a strong upflow as seen in u r . 
The maximum speed of this flow as a function of initial plasma-/3 is displayed 
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Figure 2. A series of simulations with varying plasma- (3 reveals a peak 
velocity that scales with Alfven speed. Hinode observations (scaled to these 
non-dimensional units) are consistent with these speeds. 



in Figure [2j The flows can reach speeds of up to 1/3 of the ambient Alfven 
speed, for all values of f3. Thes e values are consistent with those observed by 
Hinode/SOT ( |Berger et al^[2008l ). 

When an axial component of the field is introduced (decreasing the twist) 
the dynamics becomes richer. A series of simulations with decreasing twist 
is displayed in Figure |3j The solutions display axial Alfven waves and the 
instability transitions from a direct, growing solution to an overstable oscillation 
that eventually bifurcates to a direct mode. At smaller twists the instabilities 
are not properly captured within this 2.5D model. 



4 Summary 

We have presented a series of simulations of magnetic flux bubbles rising into 
an idealized prominence. The bubbles, driven by magnetic buoyancy, display 



morp hological characteristics similar to those observed by Hinode/SOT (Berger 
2010) and can attain rise speeds that are consistent with these observations 
The maximum speed increases in proportion to the imposed Alfven speed, even 
into the regime where the plasma- f3 is below unity. 

This work has been supported by NASA through contract NNM07AA01C. 
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Figure 3. A series of simulations in time with plasma- j3 near unity and 
increasing field twist. The background color is the azimuthal field strength, 
the lines are the axial field and the small white arrows show the flow field 
in the plane. The top panel is a purely azimuthal field ((f) = 90). The twist 
increases as <p = 89.4, 88.8 and 87.0 degrees respectively. 



